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Abstract. Two aspects of the chemical evolution of 4 He in the Galaxy 
are considered on the basis of a sample of disk planetary nebulae by the 
application of corrections due to the contamination of He from the pro- 
genitor stars. First, the He/H radial gradient is analyzed, and then, the 
helium to heavy element enrichment ratio is determined for metallicities 
up to the solar value. 



1. Introduction 

The study of radial abundance gradients and of the enrichment ratio between 
helium and the heavy elements (AY/AZ) can be performed on the basis of 
photoionized nebulae, comprising HII regions and planetary nebulae (PN, see 
for example the reviews by Peimbert 1990 and Maciel 1997). The former have 
the advantage of representing the interstellar composition directly, while for the 
latter the influence of the contamination from the progenitor stars should be 
taken into account, particularly for those elements that are dredged up to the 
outer layers of the progenitor stars, such as helium and nitrogen. 

In the case of helium, the contamination from the evolution of the pro- 
genitor stars has been recently determined on the basis of different theoretical 
models (cf. van den Hoek & Groenewegen 1997, Boothroyd & Sackmann 1999), 
so that it is now possible to apply individual corrections to the observed PN 
abundances in order to obtain improved values of the He/H radial gradient and 
of the AY/AZ ratio. In the present work, these recent calculations are used in 
order to estimate the He contamination for a sample of disk planetary nebulae 
recently studied by Maciel & Quireza (1999). The application of corrections 
to the observed abundances leads to a new determination of the He/H radial 
gradient d(Re/R)/dR and of the slope AY/AZ. 

2. The He/H radial gradient and the AY/AZ ratio 
2.1. The He/H radial gradient 

The O/H radial gradient is now well known, and can be derived for HII regions, 
planetary nebulae, hot stars, etc. It amounts roughly to dlog(0 /H) / dR ~ —0.07 
dex/kpc, and similar values have been obtained for other elements, such as 
S and Ar (Maciel 1997, Maciel & Quireza 1999). The existence of a He/H 
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gradient is much more uncertain, and results obtained both from HII regions 
and planetary nebulae show a very large uncertainty, as seen in Table 1 (see 
Esteban k Peimbert 1995 for a review). 



Table 1. Determinations of the He/H radial gradient 



Reference 


PN/HII 


d(He/H)/dfl 


dlog(He/H)/dR 






(kpc- 1 ) 


(dex/kpc) 


D'Odorico et al. (1976) 


PN 


-0.007 


-0.03 


Peimbert k Serrano (1980) 


PN 


-0.005 


-0.02 


Shaver et al. (1983) 


HII 


-0.0003 


-0.001 


Faundez-Abans k Maciel (1986) 


PN 


-0.005 


-0.02 


Pasquali k Perinotto (1993) 


PN 


-0.003 


-0.009 


Amnuel (1993) [min] 


PN 


-0.0005 


-0.002 


Amnuel (1993) [max] 


PN 


-0.006 


-0.026 


Maciel k Chiappini (1994) 


PN 


-0.0004 


-0.002 


Esteban et al. (1999) 


HII 


-0.001 


-0.004 



All determinations involving PN until now have not taken into account the 
He contamination by the progenitor stars, so that it is interesting to revise the 
values of the He/H gradient by considering the amount of He produced and 
dredged up to the outer layers of the stars and eventually shed into the nebulae. 

2.2. The AY/AZ enrichment ratio 

The helium to metals enrichment ratio AY/ AZ is generally determined adopting 
a linear variation for the helium abundance by mass with the metallicity Z of the 
form Y(Z) = Y p + (AY/AZ) Z, where Y p is the pregalactic helium abundance, as 
proposed by Peimbert and Torres- Peimbert (1974, 1976). Photoionized nebulae 
such as HII regions and blue compact galaxies are generally used (cf. Izotov 
et al. 1997, Thuan k Izotov 1998, Esteban et al. 1999), and recent work has 
also taken into account the fine structure in the main sequence of nearby stars 
(Pagel k Portinari 1998). Results are generally in the range 2 < AY/AZ < 6, 
as shown in Table 2. 

Table 2. Determinations of the AY/AZ ratio 



Reference 


object 


AY/ 


AZ 


D'Odorico et al. (1976) 


PN 


2.95 




Peimbert k Serrano (1980) 


PN 


2.2 - 


3.6 


Maciel (1988) 


PN 


3.5 ±0.3 


Pagel et al. (1992) 


HII 


6.1 ± 


:2.1 


Chiappini k Maciel (1994) 


PN/HII 


3.4- 


- 5.6 


Pagel k Portinari (1998) 


MS stars 


3±2 


Thuan k Izotov (1998) 


HII 


2.3 ± 


: 1.0 


Esteban et al. (1999) 


HII 


1.9- 


3.9 



Maciel (1988) determined Y p and AY/AZ using a sample of type II PN (Peim- 
bert 1978), and assuming that the He contamination from the central star was 
negligible. Chiappini k Maciel (1994) made a first attempt at including this 
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Figure 1. He contamination from PN progenitor stars. 

contamination in a systematic way, and added a term AY S to the ^(Z) relation, 
which can be written as 

Y = Y P + ^Z + AY S . (1) 

The stellar contribution AY S was obtained on the basis of some calculations by 
Boothroyd (private communication). The adopted values were AY S = 0.0, 0.008, 
0.015 and 0.022, which were applied to all nebulae in the sample, so that any 
differences in their individual behaviour were lost. 



3. He contamination from the PN progenitor stars 

The He excess by mass AY S can be estimated from the yields of the intermediate 
mass stars as a function of the stellar mass on the main sequence Mms an d 
total metallicity Z, adopting Z = 0.020, which is appropriate for type II PN. 
The stellar mass on the main sequence can be obtained by an initial mass-final 
mass relation as a function of the PN core mass M c , which can in principle be 
determined from the observed nebular abundances, particularly the N/O ratio. 

In this work, we have adopted the recent yields by van den Hoek & Groe- 
newegen (1997), using as comparison the calculations by Boothroyd & Sackmann 
(1999). The former includes all three dredge up processes that occur in the late 
stages of the intermediate mass stars, and generally produce larger yields than 
the latter, as can be seen in Figure 1. We have assumed that most of the He ex- 
cess that are dredged up to the surface is mixed up in the outer layers and ejected 
into the nebulae, so that our derived AY S is the largest possible correction for a 
given mass. 
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Figure 2. The adopted M c x N/O calibrations. 



In view of the uncertainties on the PN central star masses, we have adopted 
two different calibrations, referred to as high mass and low mass calibrations, 
respectively. For the high mass calibration, we adopted the M c x N/O relation 
recently proposed by Cazetta & Maciel (2000), 

M c = a + b log(N/0) + c [log(N/0)] 2 (2) 

where a = 0.689, b = 0.056 and c = 0.036 for -1.2 < log(N/0) < -0.26 and 
a = 0.825, b = 0.936 and c = 1.439 for log(N/0) > -0.26. The average initial 
mass-final mass relation was taken from the gravity distance work of Maciel & 
Cazetta (1997), and can be written as 

M c = a + oi M M s + a 2 M MS + a 3 M MS + a 4 M MS (3) 

where a = 0.5426, a x = 0.02093, a 2 = -0.01122, a 3 = 0.00447 and a 4 = 
—0.0003119. This calibration leads to core masses M c > O.67M0 and main 
sequence masses Mms > 3 M , showing a good agreement with the results from 
NLTE model atmospheres of Mendez et al. (1988, 1992). 

For the low-mass calibration we can still use equation 2, replacing the co- 
efficients by a = 0.7242, b = 0.1742 and c = for log(N/0) < -0.26. For the 
initial mass-final mass relation we have the coefficients a® = 0.4877, a\ = 0.0623, 
02 = 03 = 04 = 0. This calibration produces masses in the ranges M c > 0.55 Mq 
and Mms > 1-^© , which agree with the recent determinations of PN central 
star masses of Stasihska et al. (1997) and with the masses of type II PN origi- 
nally proposed by Peimbert (1978). Both M c x N/O calibrations are shown in 
figure 2, and the initial mass- final mass relation for the high-mass calibration 
can be seen in the figure 1 of Maciel & Cazetta (1997). 
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Figure 3. He/H abundances as a function of the galactocentric dis- 
tance for the van den Hoek and Groenewegen (1997) data with the low 
mass calibration. 



4. Results and discussion 



We have applied the corrections outlined in the previous section to a sample 
of disk planetary nebulae in order to derive the He/H radial gradient and the 
AY/AZ enrichment ratio. For details on the objects and abundances, the reader 
is referred to Maciel & Quireza (1999) and Maciel (2000). 



4.1. The He/H radial gradient 

The average He/H abundances by number of atoms are shown in Table 3, using 
the van den Hoek & Groenewegen (1997) and Boothroyd & Sackmann (1999) 
data, both for the high- and low- mass calibrations. 



Table 3. The He/H average abundances 

He/H ~ 

Uncorrected abundances 0.106 ± 0.003 

Boothroyd, low mass calibration 0.100 ± 0.003 

Boothroyd, high mass calibration 0.100 ± 0.003 
van den Hoek, low mass calibration 0.094 ± 0.003 
van den Hoek, high mass calibration 0.091 ± 0.003 



The derived gradients are negligible, and can be written as d(He/H)/di? = 
0.0000 ± 0.0004, irrespective of the calibration used, as shown for example in 
Figure 3 for the van den Hoek & Groenewegen (1997) data with the low mass 
calibration. We have adopted Rq = 7.6 kpc as in Maciel & Quireza (1999). 
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It can be seen that the correction procedure simply reduces the average 
He/H abundances, so that stars with different masses are probably scattered 
homogeneously in the whole range of galactocentric distances, thus destroying 
any systematic variations. In view of the uncertainties involved in the abun- 
dances (and distances), it is unlikely that any He/H radial gradient could be 
presently detected from planetary nebulae, and previous determinations were 
probably affected by the use of small samples. 

On the other hand, an upper limit to the He/H gradient can be derived 
on the basis of the total dispersion observed. Taking ad ~ 0.04, we have 
|d(He/H)/dR| < a d /AR ~ 0.004 kpc" 1 , or dlog(He/H)/di2 ~ -0.02 dex/kpc. 
Therefore, any existing He/H radial gradient should be lower than the O/H 
gradient by at least a factor of 3. This conclusion is in agreement with the small 
gradients derived for galactic HII regions by Esteban et al. (1999) and with 
some recent chemical evolution models for the Galaxy (Chiappini et al. 1997, 
Chiappini, this conference). 

4.2. The AY/AZ enrichment ratio 

In order to apply the correction procedure to the PN sample, we have consid- 
ered the O/H abundances as representative of the total metallicity, adopting 
Z ~ 25 O/H (Peimbert 1990, Chiappini & Maciel 1994). Since the pregalactic 
abundance can be better determined on the basis of very low metallicity objects, 
we have taken Y p as a parameter, with the values Y p = 0.23 and Y p = 0.24 (see 
for example Olive et al. 1999, Izotov et al. 1999, Steigman, this conference). 

We have taken into account all PN in our sample having metal abundances 
up to 10 6 O/H ~ 700, which corresponds approximately to the solar value, 
e(O) = log (O/H) + 12 = 8.83 (Grevesse k Sauval 1998), or Z ~ 0.017. As 
pointed out by Chiappini & Maciel (1994), He abundances of PN show some 
tendency to flatten out for very large Z, where a more sophisticated relation 
than eq. (1) would be needed. Moreover, some large metallicity PN have larger 
than average N/O ratios, so that the corrections to the He abundance are also 
larger and more uncertain. 

The main results are shown in Table 4 and in figures 4 and 5, where the 
PN are shown as filled circles. Figure 4 shows the uncorrected abundances and 
fits, and figure 5 shows the results using the corrections according to the van 
den Hoek and Groenewegen data both for the high- and low-mass calibrations. 
In both figures, the straight lines are least squares fits using Y p = 0.23 (dashed 
lines) and Y p = 0.24 (solid lines). 

For comparison purposes, figures 4 and 5 also include the Sun (©) and the 
HII regions and metal poor blue compact galaxies (empty circles) from the com- 
pilation of Chiappini & Maciel (1994). These objects have not been taken into 
account in the determination of the linear fits, and are included for comparison 
only. 

It can be seen that the scatter in the PN data is somewhat higher, which 
is partially due to the higher uncertainties in the abundances and also to the 
adopted correction procedure, so that future improvements along these lines 
would be desirable. The higher uncertainties in the PN abundances and the lower 
He/H gradient are probably responsible for the lack of a He/H x R correlation, 
while some correlation between Y and Z (or O/H) can be observed. 
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Table 4. The AY/AZ ratio 

AY/AZ 

Yj, = 0.23 



uncorrected 


5 


60 


± 





22 


Boothroyd low mass 


4 


42 


± 





20 


Boothroyd high mass 


3 


95 


± 





19 


van den Hoek low mass 


3 


59 


± 





19 


van den Hoek high mass 


2 


87 


± 





17 


Y p = 0.24 












uncorrected 


4 


73 


± 





20 


Boothroyd low mass 


3 


55 


± 





19 


Boothroyd high mass 


3 


08 


± 





17 


van den Hoek low mass 


2 


73 


± 





17 


van den Hoek high mass 


2 


01 


± 





16 




Figure 4. Uncorrected He abundances by mass Y for PN (solid cir- 
cles) as a function of O/H (top axis) and Z (bottom axis). Also shown 
are the Sun (0) and HII regions (empty circles). The straight lines are 
least squares fits for Y p = 0.23 (dashed line) and Y p = 0.24 (solid line). 
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Figure 5. The same as figure 4 using corrected abundances according 
to the van den Hoek & Groenewegen (1997) data for the high mass 
calibration (top panel) and low mass calibration (bottom panel). 
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The application of the correction procedure reduces the average He abun- 
dances and the AY/AZ ratio, and also decreases the uncertainty of the derived 
slopes. The AY/AZ ratio decreases from the range 4.7 - 5.6 to values in the 
range 2.8 - 3.6 for Y p = 0.23 and 2.0 - 2.8 for Y v = 0.24. These results are 
closer to recent independently derived ratios, as seen in Table 2 (see also Thuan, 
this conference), and to the predictions of theoretical models (Allen et al. 1998, 
Chiappini et al. 1997), particularly for the high mass calibration using the van 
den Hoek & Groenewegen (1997) data. 

The average uncertainties in the derived slopes are smaller as a consequence 
of the fact that a large number of objects has been included, with a larger 
metallicity spread. An average including both calibrations would give AY/AZ = 
3.2 ± 0.5 for Y p = 0.23 and AY/AZ = 2.4 ± 0.5 for Y p = 0.24, which can be 
compared with the uncertainties shown in Table 2. 

Finally, the apparent continuity between the low metallicity and high metal- 
licity objects of figure 5 suggests that their chemical evolution may not have been 
very different. In fact, the chemical evolution of a system is basically defined by 
its initial mass function (IMF) and star formation history (SFH). The IMF is 
now believed to be universal (Maciel & Rocha-Pinto 1998, Padoan et al. 1997). 
Blue compact galaxies have bursts of star formation, a feature that has been 
recently reinforced in our own Galaxy, as shown by Rocha-Pinto et al. (2000) 
on the basis of chromospheric ages. Therefore, the similarity in the chemical 
evolution of the different systems shown in figure 5 is probably not surprising. 
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